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Abstract 



A unit disk graph is the intersection graph of n congruent disks 
in the plane. Dominating sets in unit disk graphs are widely stud- 
ied due to their application in wireless ad-hoc networks. Because the 
minimum dominating set problem for unit disk graphs is NP-hard, 
£^ \ numerous approximation algorithms have been proposed in the liter- 

C/5 ■ ature, including some PTAS. However, since the proposal of a linear- 

i O i \ time 5-approximation algorithm in 1995, the lack of efficient algorithms 

attaining better approximation factors has aroused attention. We in- 
£N| \ troduce a linear-time 0(n + m) approximation algorithm that takes 

the usual adjacency representation of the graph as input and outputs 
a 44/9- approximation. This approximation factor is also attained by 
a second algorithm, which takes the geometric representation of the 



\ graph as input and runs in 0(n\ogn) time regardless of the number 



of edges. Additionally, we propose a 43/9-approximation which can 
be obtained in 0(n 2 m) time given only the graph's adjacency repre- 
. sentation. It is noteworthy that the dominating sets obtained by our 

algorithms are also independent sets. 



^ ! 1 Introduction 



A unit disk graph G is a graph whose n vertices can be mapped to points in 
the plane and whose m edges are defined by pairs of points within Euclidean 
distance at most 1 from one another. Alternatively, one can regard the 
vertices of G as mapped to coplanar disks of unit diameter, so that two 
vertices are adjacent whenever the corresponding disks intersect. 
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A dominating set D is a subset of the vertices of a graph such that 
every vertex not in D is adjacent to some vertex in D. An independent 
dominating set is a dominating set which is also an independent set. Note 
that any maximal independent set is an independent dominating set. 

Dominating sets in unit disk graphs are widely studied due to their ap- 
plication in wireless ad-hoc networks |12j . Since it is NP-hard to compute 
the minimum dominating set of a unit disk graph [2] , several approximation 
algorithms have been proposed G3 [TOl 021 [151 GO]. Such algorithms 
are of two main types. Graph-based algorithms receive as input the adja- 
cency representation of the graph and assume no knowledge of the point 
coordinates, whereas geometric algorithms work in the Real RAM model of 
computation and receive solely the vertex coordinates as inpulQ. Thus far 
these two types of algorithms have been tackled separately in the literature 
for the dominating set problem in unit disk graphs. In this paper, we in- 
troduce approximation algorithms of both types, benefiting from the same 
approximation factor analysis. 

Previous algorithms A graph-based 5-approximation algorithm that runs 
in 0(n + m) time was presented in [12j . The algorithm computes a maximal 
independent set, which turns out to be a 5-approximation because unit disk 
graphs contain no as induced subgraphs. 

Polynomial-time approximation schemes (PTAS) were first presented as 
geometric algorithms [10J and later as graph-based algorithms [15]. Also, 
a graph-based PTAS for the more general disk graphs was proposed in [9]. 
Unfortunately, the complexities of the existing PTAS are high degree poly- 
nomials. For example, the PTAS presented in [15] takes 0(n 225 ) time to 
obtain a 5-approximation (using the analysis from [3j). Although its anal- 
ysis is not tight, the running time is too high even for moderately large 
graphs. The reason is that these PTAS invoke a subroutine that verifies (by 
brute force) whether a graph admits a dominating set with k vertices. Such 
subroutine is applied to several subgraphs, and the value of k grows as the 
approximation error decreases. A similar strategy was used in [llj to obtain 
a PTAS for the minimum independent dominating set. 

The lack of fast algorithms with approximation factor less than 5 was re- 
cently noticed in [3] , where geometric algorithms with approximation factors 
of 4 and 3 and running times respectively 0(n 9 ) and 0(n 18 ) were presented. 
While a significant step towards approximating large instances, those algo- 
rithms require the geometric representation of the graph, and their running 
times are still polynomials of rather high degrees. Linear and near-linear- 
time approximation algorithms constitute an active topic of research, even 
for problems that can be solved exactly in polynomial time, such as maxi- 

1 The Real RAM model is a technical necessity, otherwise storing the coordinates of the 
vertices would require an exponential number of bits [13] . 
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mum flow and maximum matching [IJ [19] . 

It is useful to contrast the minimum dominating set problem with the 
maximum independent set problem. While a maximal independent set is a 
5-approximation to both problems, it is easy to obtain a geometric 3-approx- 
imation to the maximum independent set problem in 0(n log n) time |14j . 
In the graph-based version, a related strategy takes roughly 0(n 5 ) time, 
though. No similar results are known for the minimum dominating set 
problem. 

Packing problems usually arise in situations where one wants to enclose 
non-overlapping objects into a recipient of given shape covering the greatest 
possible fraction of the recipient area. Unit disk graphs are subject to a 
number of packing constraints that limit the size of independent sets (which 
correspond to disjoint disks) as a function of the distance between vertices. 
The existing PTAS for dominating sets in unit disk graphs are based on some 
of these packing constraints, such as the bounded growth property: the size 
of an independent set formed by vertices within distance r of a given vertex 
is at most (1 + 2r) 2 . Note, however, that the bounded growth property is 
not tight. For example, for r = 1, it gives an upper bound of 9 vertices 
where the actual maximum size is 5. Since the bounded growth property is 
strongly connected to the problem of packing circles in a circle [7], obtaining 
exact values for all r seems unlikely. 

Our contribution Our main result consists of the two approximation 
algorithms given in Section a graph-based algorithm, which runs in linear 
0(n + m) time, and its geometric counterpart, which runs in 0(n log n) time 
in the Real RAM model, regardless of the number of edges. 

The approximation factor of both algorithms is 44/9, and the strategy 
in both cases is to construct a 5-approximate solution using the algorithm 
from [12], then perform subsequent local improvements to that initial dom- 
inating set. Our main lemma (Lemma [9|) uses forbidden subgraphs to show 
that a solution that admits no local improvement is a 44/9-approximation. 
Since the dominating sets produced by our algorithms are independent sets, 
the same approximation factor holds for the independent dominating set 
problem. 

Proving that a certain graph is not a unit disk graph (and is therefore a 
forbidden induced subgraph) is no easy featd. We make use of an assortment 
of results from discrete geometry in order to prove properties of unit disk 
graphs that are interesting per se. For example, we use universal covers and 
disk packings to show that the neighborhood of a clique in a unit disk graph 
contains at most 12 independent vertices. These properties, along with a 
tighter version of the bounded growth property, are collected in Section [2J 

2 The fastest known algorithm to decide whether a given graph is a unit disk graph is 
doubly exponential [17] . 
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and allow us to show that certain graphs are not unit disk graphs. Conse- 
quently, the analyses of our algorithms employ a broader set of forbidden 
subgraphs, including, but not being limited to, the Kifi. 

Additionally, in Section HI we show that a possible, somewhat natu- 
ral refinement to our graph-based algorithm leads indeed to a tigher 43/9- 
approximation, albeit for the price of an extra 0(n 2 ) multiplying factor in 
its time complexity. 

2 Forbidden subgraphs 

In this section, we introduce some lemmas about the structure of unit disk 
graphs. These lemmas will be applied to prove our approximation factors in 
Sections [3] and |H We start by stating three previous results from the area 
of discrete geometry. The first lemma comes from the study of universal 
covers (for a recent survey see [8]). 

Lemma 1 (Pal |16j). If a set of points P has diameter 1, then P can be 
enclosed by a circle of radius l/\/3- 

Packing congruent disks in a circle is a well studied problem. Exact 
bounds on the radius of the smallest circle packing k congruent disks are 
known for some small values of k, namely k < 13 and k = 19 [7]. The bound 
for k = 13 will be useful to us. 

Lemma 2 (Fodor [7j). The radius of the smallest circle enclosing 13 points 
with mutual distances at least 1 is (1 + y/E)/2. 

The density of a packing is the ratio between the covered area and the 
total area. The following general upper bound is useful when no exact bound 
is known. 

Lemma 3 (Fejes Toth [6]). Every packing of two or more congruent disks 
in a convex region has density at most ir/y/12. 

Given a graph G = (V, E) and a vertex v S V, let N(v) denote the open 
neighborhood of v and let N[v] = N(v) U {v} denote the closed neighborhood 
of v. More generally, the open r -neighborhood of a vertex v is the set of 
vertices w such that the distance between v and w in G is exactly r, while 
the closed r -neighborhood of a vertex v is the set of vertices w such that 
the distance between v and w in G is at most r. For a set S C V, we let 
Ns(v) = N(v) n S and N s [v] = N[v] n S. Finally, given a subgraph H of 
G, the closed neighborhood of H, denoted N[H], is the set of vertices that 
belong to the closed neighborhood of some vertex of H. The following two 
lemmas concern neighborhoods in unit disk graphs. 

Lemma 4. The closed neighborhood of a clique in a unit disk graph contains 
at most 12 independent vertices. 
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Proof. By LemmaHJ the points which define a clique in a unit disk graph are 
contained inside a circle of radius 1/V3. Therefore, the points corresponding 
to the closed neighborhood of such clique are contained inside a circle of 
radius 1 + (l/y/3). By Lemma [21 we have that a circle enclosing 13 points 
with mutual distances at least 1 has radius at least (1 + y/5)/2. Since 
(l + V5)/2 > l + (l/\/3), the lemma follows. □ □ 

Lemma 5. Given an integer r > 1, the closed r -neighborhood of a vertex in 
a unit disk graph contains at most [ir(2r + l) 2 /\/l2| independent vertices. 

Proof. All n disks of diameter 1 corresponding to the closed r-neighborhood 
of a vertex v must be enclosed by a circle Z of radius (2r + l)/2 centered 
on v. Each disk of diameter 1 has area tt/4 and Z has area (2r + l) 2 7r/4. 
Using Lemma[3j we have (n 7r/4)/((2r + l) 2 7r/4) < 7r/\/l2, and the lemma 
follows. □ □ 

We say that a graph G is (k,£)-pendant if there is a vertex v in G with 
k vertices of degree 1 in the open neighborhood of v and t vertices of degree 1 
in the open 2-neighborhood of v. We refer to v as a generator of the (k, £)- 
pendant graph. The following lemma bounds the value of the parameter £ 
for (4, £)-pendant unit disk graphs. 

Lemma 6. If G is a (4,£)-pendant unit disk graph, then t < 8. 

Proof. Let v be a generator of G. Since is a forbidden induced sub- 
graph [12] and v has 4 neighbors of degree 1, we have that the remaining 
neighbors of v together with v itself form a clique. By Lemma 0J it follows 
that 4 + I < 12. □ □ 

Next, we consider the case of (3, £)-pendant unit disk graphs. 

Lemma 7. If G is a {?>,£) -pendant unit disk graph, then £ < 16. 

Proof. Let v be a generator of G. Since two vertices are adjacent if and only 
if their Euclidean distance is at most 1, the closed neighborhood of v lies 
inside a circle of radius 1 centered at v. Let u be a neighbor of v with 
degree 1. We divide the aforementioned circle into 6 congruent sectors 
si, . . . , sq in such a way that u sits on the boundary of two adjacent sectors 
s\,S2- Since the diameter of each sector is 1 and u has degree 1, we have 
that s\ U S2 only contains the single vertex u. By the same argument, the 
remaining two neighbors of v that have degree 1 are contained in 2 of the 
remaining 4 sectors, which we call ss,s^ (note that the sectors do not nec- 
essarily appear in the order s±, . . . , se). Therefore, only sectors S5, sq may 
contain the neighbors of v that have degree greater than 1. 

Notice that the neighbors of the vertices in S5 U sq should be located 
within Euclidean distance at most 1 from 55 U sg. Consequently, the disks of 
unit diameter corresponding to such vertices should be completely contained 
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Figure 1: Upper bounds to the area of the Minkowski sum of S5 U sq and a 
disk of radius 1.5 in two different scenarios. 

inside the region defined as the Minkowski sum of S5 U sg and a disk of radius 
1.5, i.e. the region within distance at most 1.5 from 55 U sq. We compute 
upper bounds to such area by considering two cases, depending on whether 
55 and sq are opposite sectors (bounded by the same pair of straight lines) 
or not. 

First, we consider the case when 55 and sq are not opposite sectors. In 
this scenario, S5 U sq is contained in a semicircle of radius 1 centered at the 
generator v, as represented in Figure [TJa). We define a region R as the locus 
of points within distance at most 1.5 from the aforementioned semicircle of 
radius 1. The area a of R is therefore a = 3 + 177r/4. 

By Lemma El the number of disks of unit diameter contained in R is at 
most a ■ (y^) ' (774) < 18 - 8814 < 19. Since k = 3 of these disks are the 
degree- 1 neighbors of v, we have i < 16. 

Second, we consider the case when S5 and sq are opposite sectors. In this 
scenario, the region S5 U sq is not convex. In order to obtain a convex region 
R' , we define R' as the locus of points within distance at most 1.5 from 
the convex hull of S5 U sq (see Figure QJb)). The area a' of R' is therefore 
7^3/2 + 43vr/12. 

By Lemma [3l the number of disks of unit diameter contained in R' is at 
most a' ■ (jfi) ■ [yjij < 19.9989 < 20. Since k = 3 of these disks are the 
degree- 1 neighbors of v, we have i < 17. □ □ 

Finally, the following lemma holds for the general case. 

Lemma 8. If G is a (k, I) -pendant unit disk graph, then k + 1 < 22. 

Proof. Immediately from Lemma [5] with r = 2. □ □ 

3 Linear-time 44/9-approximation 

In this section, we present two 44/9-approximation algorithms. The key 
property to analyze the approximation factor is presented in Lemma[9l while 
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the running time analyses are presented in Sections 13,11 and 13,21 

Hereafter, let G = (V, E) be a unit disk graph, and let D C V be an 
independent dominating set of G. If v G D and uv G E, we say that v 
dominates u and, conversely, that u is dominated by v. 

As already mentioned, unit disk graphs are free of induced K\£. There- 
fore, at most 5 vertices of D may belong to the closed neighborhood of any 
given vertex v G V. A corona is a set C C D consisting of exactly 5 neigh- 
bors of some vertex c G V \ D. Such a vertex c, which is not necessarily 
unique, is called a core of the corona C, whereas the 5 vertices of the corona 
are referred to as the corona's petals. Notice that the subgraph induced by 
a corona C and a corresponding core c is a star, i.e. a graph formed by an 
independent set and a universal vertex. 

A corona C is said to be reducible if there is a core c of C such that 
DL){c}\C is a dominating set. If no such core exists, C is dubbed irreducible. 
Given a reducible corona C and a corresponding core c, we refer to the 
operation that converts D into the smaller dominating set D U {c} \ C as a 
reduction of C with respect to c. 

Lemma 9. Let G = (V, E) be a unit disk graph, D an independent domi- 
nating set in G, and D* a minimum dominating set of G. If all coronas in 
D are irreducible, then p = \D\/\D*\ < 44/9. 

Proof. We use a charging argument to bound the ratio between the cardi- 
nalities of D and D* . Consider that each vertex u G D splits a unit charge 
evenly among the vertices in the closed neighborhood Nd* [u]. The function 
/ : D* — ^ (0, 5] below corresponds to the total charges assigned to each ver- 
tex v* G D* , accumulating the (fractional) charges that v* received from the 
vertices in Nr)[v*]: 



ueN D [v*} 1 u L Jl 

Note that, since D and D* are dominating sets, neither Nd*[u] nor Nd[v*] 
are ever empty, and f(v*) < Nd[v*]. Such function / allows us to write the 
cardinality of D as 

i^i = E 

v*eD* 

Since 

1^1 £„. e p*/fr*) 

is precisely the average value of /(•) over the elements of D* , we obtain 
the desired bound p < 44/9 by showing that the existence of vertices c* 
in D* with /(c*) > 44/9 is counterbalanced by a sufficiently large number 
of vertices r* in D* with f(r*) < 4. 
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Figure 2: Figure for the proof of Lemma [H A proper, induced subgraph, 
where squares were used for a subset of D* , solid circles for the corona 
CCD, and hollow circles for vertices not in D U D*. Vertices w and r* are 
respectively witness and reliever of core c*. 

Before we continue, note that f(c*) > 44/9 means exactly f(c*) = 5, 
because the summation in ([1]) has at most 5 terms, all of which are of the 
form 1/i for integer i > 1. 

Thus, let c* be a vertex in D* with /(c*) = 5. Clearly, c* ^ D, otherwise 
/(c*) < | iVo [c* ] I = 1, because D is an independent set. Moreover, c* must 
have exactly 5 neighbors in D, since a greater number of neighbors in D 
would imply the existence of an induced K\ § in G, which is not possible, 
and a lesser number would imply /(c*) < |7Vd[c*]| < 4, a contradiction. 
Vertex c* is therefore a core. 

Now let C C D be the corona of which c* is a core. Since C is irreducible 
(by the hypothesis of the lemma), there must be a vertex w £ V\(CL){c*}), 
such that: 

(i) w is not adjacent to c*; 

(ii) w is only dominated, in D, by vertices that belong to C. 

We call w a witness of c*, meaning the corona having c* as a core cannot 
be reduced with respect to c* due to the existence of w. Note thatR' 

(iii) for all u G C, the only vertex in iVD*[it] is the very core c* of C, 

otherwise the contribution of some u £ C in the summation yielding f(c*) — 
given by flU — would be at most 1/2, and /(c*) would be at most 9/2 < 5, 
a contradiction. In particular, fact (iii) above implies that the witness w, 
which is adjacent to at least one vertex in C, cannot belong to D*. But D* 
is a dominating set, so there must exist a vertex r* G D* that is adjacent 
to w, and r* ^ c* due to (i). We call r* a reliever of c*. Figure [2] illustrates 
this situation. 

We now show that |iYD[r*]| < 4. For the sake of contradiction, assume 
| JVd [r*] | > 4. Because G is free of induced K% t Q, such number must be 
exactly 5, so that r* is the core of a corona C C D. Such corona must be 
disjoint from corona C, otherwise there would be a vertex in CD C adjacent 
to more than one vertex in D* , namely c* and r* , contradicting (iii). Now, 
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due to (ii), we have Nc(w) = 0, hence C U {w} is an independent set 
of G, constituting, along with the core r* , an induced K\$ in G. This is a 
contradiction, because G is a unit disk graph. Thus, |A^[r*]| < 4. Since 
/(r*) < |iVo[r*]|, we have f(r*) < 4. 

We have just shown that the existence of a vertex c* in D* with f(c*) = 5 
implies the existence of a vertex r* £ D* such that /(r*) < 4. Were 
this correspondence one-to-one, we would be able to state that the average 
of /(•) over the elements of D* was no greater than 4.5. Unfortunately, this 
correspondence is not necessarily one-to-one, as exemplified by the graph 
in Figure El for which a disk model is given in Figure 0] with coordinates 
presented in Table [TJ 

Still, the lemmas in Section [2] allow us to bound the ratio between the 
number of vertices c* with f(c*) = 5 and the number of vertices r* for which 
the values of / are significantly lower. Let r* € D*\D be a reliever. In order 
to obtain the claimed bound, we consider two cases according to the size 
of N D [r*]: 

(i) I-^YdI?"*] I < 3. By Lemma[5l the closed 4-neighborhood of r* contains 
at most 73 independent vertices. Since each corona contains 5 independent 
vertices, at most [(73 — |-/Vd[t*]|)/5J = 14 cores may share a common re- 
lieve^. Let c*,...,c* 4 G D* denote such cores. If l-^Vofr*]) < 3, then the 
average value of /(•) among r*, d[, . . . , c* 4 is at most 

1-3 + 14-5 

-— < 4.867. 

15 

(ii) |AT D [r*]| = 4. By Lemma El if |A?i)[r*]| = 4, then at most 8 cores 
c*,...,Cg may have r* as their common reliever, for otherwise we obtain 
a (4, 9)-pendant graph, which cannot be a unit disk graph. Thus, the average 
value of /(•) among r*, c|, . . . , c| is at most 

1-4 + 85 = 44/9 = 4.888 . . . 

The worst case is therefore the one in which |-/V£>[r*]| = 4, for an average 
p = 44/9, and the lemma follows. □ □ 

3.1 Graph-based algorithm 

By Lemma [9l an independent dominating set with no reducible coronas is a 
44/9- approximation to the minimum dominating set. We now describe how 
to obtain such set in linear time given the adjacency list representation of 
the graph. 

We can easily compute a maximal independent set D, which is a 5- 
approximation to the minimum dominating set [12], in 0(n + m) time. An 

3 We would like to thank an anonymous referee for this simplified argument. 
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Figure 3: Unit disk graph where 4 distinct cores c*, . . . , c| share the same 
reliever r*. 



w 2 
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Figure 4: Geometric model for the graph in Figure O Due to scale/precision 
issues, some disks appear to touch one another when in fact they do not. For 
clarity, disks that actually touch one another have their centers connected 
by a straight line. Moreover, there are no concentric disks among those for 
the witnesses w%, . . . ,w^. The coordinates of the centers, comprising only 
integer numbers, are given in Table [TJ 
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r* 


(0,0); 




r* 


(-2492384,879081), 


W\ 


r* 


(-1310377,2686162), 


W 2 


r* 


(1310377, 2686162), 


W 3 


n* 


(2492384,879081), 


W\ 



-492423,870355), 
-484809,874619), 



remaining vertices: 
(±776025,3531423), 
(±2309705,2722805), 
(±1776763,3570742), 
(±503019,-864274), 



(±1492384,879081), 
(±3491646,917468), 
(±1840296,1838114), 
(±2957226,1764474), 



(±999986,5235), 
(±3023782,31960), 
(±2022913,-3866), 
(±810377, 1820137). 



Table 1: Coordinates of the centers of the disks in Figure [U where all 
diameters are equal to 1000001. 



independent dominating set with no reducible coronas can then be obtained 
by iteratively performing reductions. However, naively performing such re- 
ductions leads to a running time of 0(n 2 m), since (i) there are 0(n) candi- 
dates to being the core of a reducible corona, (ii) detecting whether a vertex 
v is in fact the core of a reducible corona by inspecting the 3-neighborhood 
of v takes 0(m) time, and (hi) we may need to reduce a total of 0(n) coro- 
nas. Fortunately, the following algorithm modifies the set D and returns an 
independent dominating set with no reducible coronas in 0(n + m) time. 

(1) For each vertex v G V \ D, compute Nd(v). 

(2) For each vertex v £ V \ D, if \N D (v)\ = 5, add N D (v) to the list of 
coronas C (unless it is already there). 

(3) Let B <— 0. For each corona C G C, if there is a vertex c such that 
D U {c} \ C is a dominating set, then add c to the set B. 

(4) Choose a maximal subset B 1 of B such that the pairwise distance of the 
vertices in B' is at least 5. 

(5) For each vertex c 6 B' , perform a reduction D <— DU {c} \ Nu(c). 

(6) Repeat all the steps above until B' = 0. 

The algorithm is correct since all changes made to D along its execution 
preserve the property that D is an independent dominating set. Notice that, 
in step (4), we only reduce coronas that are sufficiently far from each other, 
in order to guarantee that we do not reduce a corona that may have ceased 
to be reducible due to a previous reduction. Moreover, the algorithm always 
terminates because the size of D decreases at every iteration, except for the 
last one. Next, we show that the running time is 0(n ± m). 
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Step (1) can be easily implemented to run in 0(n + m) time. To execute 
step (2) in 0(n + m) time, we must determine in constant time whether a 
corona is already in the list C. This can be achieved by indexing each corona 
C by an arbitrary vertex v G C (say, the one with the lowest index), and 
by storing with v a list of coronas that are in C and whose index is v. Note 
that, because of the packing constraints inherent to unit disk graphs, the 
number of coronas that contain a given vertex is 0(1). 

Step (3) can be implemented as follows (for each corona C G C): 

(3a) Let S\ be the union of the open neighborhoods of the 5 petals of C. 

(3b) Let 52 be the subset of all vertices w G Si with Njj{w) C C . 

(3c) Let 53 be the intersection of the closed neighborhoods N[v] of all v G 
S 2 U C. 

(3d) If S3 7^ 0, then add an arbitrary vertex of S3 to the set B. 

The steps above take 0(n + m) total time when executed for all coronas 
C G C, because the number of coronas that contain or are adjacent to a 
given vertex is also O(l) by packing constraints. 

It is easy to perform steps (4) and (5) in linear time. It remains to show 
that the whole process is only repeated for a constant number of iterations. 
Let Bi, . . . , Bj. denote the set of reducible coronas at each iteration of the 
algorithm with = 0. Note that the reductions performed in step (5) never 
create a new reducible corona. Therefore B\ D • • • D B^. Let C denote a 
corona that was reduced in the last iteration k. If C was not reduced during 
a previous iteration i < k, then another corona within distance 5 from C 
was reduced at that very iteration i. Since, again by packing constraints, 
the maximum number of coronas within constant distance from C is itself a 
constant, we have k = 0(1). 

The following theorem summarizes the result from Section 13.11 

Theorem 10. Given the adjacency list representation of a unit disk graph 
with n vertices and m edges, we can find a 44/ '9 -approximation to the min- 
imum dominating set in 0(n + m) time. 

3.2 Geometric algorithm 

In this section, we describe how to obtain an independent dominating set 
with no reducible corona in 0(n log n) time given the geometric representa- 
tion of the graph. The input is therefore a set P of n points. Without loss of 
generality, we assume that the corresponding unit disk graph is connected 
(otherwise, we can compute the connected components in 0(n log n) time 
using a Delaunay triangulation [3]). We use terms related to vertices of the 
graph and to the corresponding points interchangeably. For example, we 
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say a set of points is independent if all pairwise distances are greater than 
1. 

We want the points of P to be structured in suitable fashion. Thus, as a 
preliminary step, we sort the points by x-coordinates and by y-coordinates 
separately (such orderings will also be useful later on), and we partition 
the points of P according to an infinite grid with unitary square cells by 
performing two sweeps on the sorted points. Without loss of generality, we 
assume that no point lies on the boundary of a grid cell. Given p G P, let 
cr(p) denote the grid cell that contains p. We refer to the set of at most 8 
non-empty grid cells surrounding a cell Q as the open vicinity of Q, denoted 
N(Q), and to the union of Q and its open vicinity as the closed vicinity of 
Q, denoted N[Q]. Note that a point p can only be adjacent to points in the 
closed vicinity of o~(p), that is, N\p] C N[a(p)]. Each point p G P stores a 
pointer to its containing cell o~(p). Also, each cell stores the list of points it 
contains and pointers to the cells in its open vicinity. Since the diameter of 
the point set is at most n due to the graph connectivity, this whole step can 
be done in 0(n log n) time. 

We are now able to show how to compute a maximal independent set D 
efficiently We begin by making a copy P' of P, and by letting D <— 0. Then 
we repeat the two following steps while set P' is non-empty, (i) Choose an 
arbitrary point p G P' and add it to set D. (ii) For each point jf in the 
closed vicinity of o~(p), remove p' from P' if \\pp'\\ < 1. When P' becomes 
empty, D is an independent dominating set. This process takes 0(n) time 
due to the two following facts. First, a cell belongs to the closed vicinity of 
a constant number of cells. Second, the maximum number of points inside 
a cell with pairwise distances greater than 1 is also a constant. 

We now have that D is a maximal independent set, and therefore 
a 5- approximation to the minimum dominating set. Next, we show how 
to modify D in order to produce an independent dominating set with no re- 
ducible corona, therefore a 44/9- approximation to the minimum dominating 
set. The algorithm follows a close parallel to the one in Section T3. 11 but each 
step takes no more than 0(n log n) time using the geometric representation 
of the graph. 

Since D is an independent set and a grid cell Q has side 1, a simple 
packing argument shows that \D n Q\ < 4. We store the set D n Q in the 
corresponding cell Q. In order to compute Nu(p), it suffices to inspect at 
most the 36 points in D n Q for Q G N[o~(p)]. We can then build a list of 
coronas in 0(n) time (steps (1) and (2) of Section IBTTj) . 

To perform step (3), we need to find out whether there is a core c such 
that DU{c}\C is a dominating set, for each corona C = {pi, . . . ,ps}- First, 
we make Si the union of Nz>(pi) for 1 < i < 5. Then, we let 52 be the subset 
of Si containing only the points p with No{p) C C. These first two steps 
are similar to steps (3a) and (3b) in Section 13.11 The remaining sub-steps 
of step (3) are significantly different, though. 
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We proceed by making £3 = S2UC. We need to determine whether there 
is a point p £ S3 that is adjacent to all points in S3. For each p £ S3, let j3{p) 
denote the disk of radius 1 centered at p. Let R denote the convex region 
defined by the intersection of f5(p) for all p € H. A point p is adjacent 
to all points in S3 if and only if p £ R. We can compute the region R 
in O ( I S*3 1 log I S3 1) time using divide-and-conquer in a manner analogous to 
half-plane intersection [3]. We can then test whether each point p £ S3 
belongs to the region R in logarithmic time using binary search (remember 
the points were previously sorted). If there is at least one point p £ S3 Pi R, 
then we add p to the set B. Therefore, the whole step (3) takes 0(n log n) 
time. 

In step (4) of the geometric algorithm, we choose an alternative set 
B' C B which can be computed in 0(n) time as follows. For each p £ B, we 
add p to B' and then remove from B all points that are contained in the cells 
within Euclidean distance at most 4 of cr(p). Since by packing constraints 
there are O(l) points in the intersection of D and the closed vicinity of a 
cell, we can easily perform step (5) in 0(n) time. 

The following theorem summarizes the result from Section 13.21 

Theorem 11. Given a set of n points representing a unit disk graph, we 
can find a 44/ '9 -approximation to the minimum dominating set in 0(n log n) 
time in the Real RAM model of computation. 

4 Achieving a 43/9-approximation 

In the previous section, a 44/9- approximation was obtained by reducing 
coronas of a maximal independent set D of graph G, that is, by subsequently 
replacing 5 petals with 1 core in D as long as that operation preserved 
dominance. A natural step to tighten the approximation factor is to allow 
for weak reductions, whereby the 5 petals of a corona C are removed from 
the independent dominating set D, yet not only a core c of C is inserted 
into D but also some mutually non- adjacent witnesses of c, as long as their 
number is no greater than 3 and they dominate all witnesses of c. By doing 
so, the weak reduction of C (with respect to c) preserves dominance and still 
shaves off at least one unit from the size of D. If such operation is possible 
on a corona C, then C is said to be weakly reducible. A core c which has 
4 (or more) mutually non- adjacent witnesses is said to be an overwhelmed 
coreJll 

4 If c is an overwhelmed core of a corona C, then it might be the case that a weak 
reduction on C with respect to c is still possible. If the subgraph G[W] induced by the set 
W of witnesses of c admits an independent dominating set W' C W of size no greater than 
3, then (D \ C) U {c} U W' is still an independent dominating set of G, and its cardinality 
is strictly less than \D\. However, one cannot decide efficiently whether such dominating 
set exists, and that will not be required by our algorithm. 
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We consider the graph-based algorithm presented in Section 13,11 with 
some modifications to cope with weak reductions. The whole modified al- 
gorithm becomes: 

(1) For each vertex v G V \ D, compute Njj(v). 

(2) For each vertex v G V \ D, if \N D (v)\ = 5, add C = N D (v) to the list of 
coronas C (unless it is already there), and add v to the list Lc containing 
the cores of C. 

(3) Let B be an initially empty mapping of cores onto sets of witnesses. 
For each corona C G C, if there is a vertex c G Lc and an independent 
set W c C V with at most 3 vertices such that (D \C) U {c} U W c is a 
dominating set, then add (c, W c ) to B. 

(4) Choose a maximal subset B 1 of the cores in B such that the pairwise 
distance of the vertices in B' is at least 5. 

(5) For each vertex c G B', perform a (weak) reduction D <— (D \ Nd(c)) U 
{c} U W c . 

(6) Repeat all the steps above until B' = 0. 

The new step (3) can be implemented as follows (for each corona C G C): 

(3a) Let Si be the union of the open neighborhoods of the 5 petals of C. 

(3b) Let 5*2 be the subset of all vertices w G S\ with Nd(w) C C. 

(3c) For each core c G Lc, greedily obtain a maximal independent set Wc 
of S 2 \ N[c}. If \W C \ < 3, add (c, W c ) to £ and break. 

Because each core is now evaluated separately as to whether or not it al- 
lows for a weak reduction (whereas in the algorithm of Section ^. ll a constant 
number of set operations per corona was executed), step (3c) now dominates 
the complexity of the whole sequence of steps (1) to (5), with a total 0(nm) 
time for running on all coronas C G C. Another important difference, as 
far as time complexity goes, is that, in this modified algorithm, the num- 
ber of iterations of the main loop — steps (1) to (5) — is no longer O(l), 
due to the fact that new (weakly) reducible coronas can be created, as illus- 
trated in FigureO However, such number is certainly 0(n), because the size 
of D decreases by at least one unit at each iteration. Hence the overall time 
complexity of the modified algorithm is 0(n 2 m). 

Next, we establish the approximation factor of the modified algorithm. 
Note that step (3c) asserts that the coronas that are not (weakly) reduced 
by the algorithm present only overwhelmed cores. 
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Figure 5: Example of graph with a maximal independent set (represented 
by solid circles) containing only 1 weakly reducible corona. After every weak 
reduction performed by the modified algorithm, except for the last one, a 
new weakly reducible corona is created. 

Lemma 12. Let G = (V, E) be a unit disk graph, D an independent dom- 
inating set in G, and D* a minimum dominating set of G. If all coronas 
in D have only overwhelmed cores, then p = \D\/\D*\ < 43/9. 

Proof. We follow a strategy similar to that of Lemma [9j by employing the 
same function / : D* — >■ (0, 5] defined in (pQ), we prove there is an appropriate 
balance among vertices with high (> 43/9) and low (< 4) images under /, 
thus yielding an average value for /(•) that is no greater than 43/9 — the 
claimed approximation factor. 

Let c* be a vertex in D* with /(c*) > 43/9. Because 43/9 > 4.5, and all 
(at most 5) terms in the summation that yields / are either 1 or no greater 
than 0.5, we have that f(c*) > 43/9 implies f(c*) = 5, so c* is a core. By 
hypothesis, c* is overwhelmed, hence c* possesses at least 4 mutually non- 
adjacent witnesses w\, . . . , W4, each one implying the existence of a reliever 
r* £ D* with f{r*) < 4. Surely, such relievers need not be all distinct, 
and different cores may share a common reliever. Still we can define a 
suitable core-to-reliever relationship and determine an upper bound to the 
average value of /(•) in D* — therefore to the approximation factor — by 
investigating the cardinalities of the sets in such relationship. 

Suppose there are t cores c* £ D* such that f(c*) = 5. For i = 1, . . . , t, 
we let Ci be the corona having c* as a core, Wi be the set of (at least 4) mu- 
tually non-adjacent witnesses of c*, and Ri be the set of relievers of c*. Now, 
we construct a bipartite multigraph G' = (C U R 1 , E') as follows. The parts 
of G' are C' = {c*, ...,c£} and R' = R\ U ... U Rt- The multiset E' contains, 
between each core c* £ C' and reliever r| £ R', a number (f)(c*, rj) of parallel 
edges that is equal to the number of petals of Ci adjacent to witnesses w 
(of c*) such that w is a neighbor of r*. Thus, the total number of edges 
incident to a core c* £ C is 

d(c*)= <K<$> r i)> 
Analogously, the total number of edges incident to a reliever r* £ R' is 

c*ec> 
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We now obtain an upper bound p' = 43/9 for the average value of /(•) 
over C U R' . Observe that C contains all vertices c* of D* such 
that f(c*) > 43/9, hence the average value p of /(•) over the whole set 
D* 5 C U R' cannot be any greater. 

Of course the average value we are interested in depends on the core-to- 
reliever ratio |C"|/|i?'| in G': the more cores (respectively, the less relievers) 
in C (in R'), the greater the average. Therefore, in order to obtain the 
desired upper bound p' , we must consider the case in which the elements 
of C (respectively, of R') have degrees in G' that low (as high) as 

possible. 

First, notice that, if c* is a core of corona C and /(c*) = 5, then it is 
not possible that more than 2 witnesses of c* sharing a common reliever r* 
are adjacent to the same petal p S C. Otherwise, let w\,W2 and W3 be such 
witnesses. Since p and r* are non-adjacent (due to the image of c* under / 
being 5), the subgraph of G induced by {p, r*, w%, W2, W3} is a i^2,3- This is 
a contradiction, because the i^2,3 is not a unit disk graph [IS] . Thus, since 
c* has at least 4 witnesses, the number of petals of C adjacent to witnesses 
of c* — and therefore the degree of each core in C — is at least [4/2] =2. 

Such lower bound can be improved for cores c* G C having a reliever r* 
with exactly 4 neighbors in D. If a reliever r* has 4 neighbors in (the inde- 
pendent set) D, then the remaining neighbors of r* form a clique. Now, if a 
witness w of c* is one such neighbor of r*, then the other (at least) 3 mutually 
non-adjacent witnesses of c* must have relievers distinct from r* . Moreover, 
those witnesses must be adjacent to at least [3/2] = 2 petals of the corona 
having c* as a core. This means there are at least 2 edges connecting c* to 
its relievers in R' \ {r*}, plus 1 edge connecting c* to r* , hence the degree 
of c* in G' is at least 3. 

As for the maximum degree of vertices r* £ R' , they basically depend 
on the number Nu(r*) of neighbors of r* in D, and now we employ some 
of the geometric lemmas of Section [2] to infer suitable upper bounds. Recall 
that any set of petals is an independent set. 

• If Nu(r*) = 4 (implying f(r*) < 4), then, by Lemma [6] there are at 
most 8 petals in the 2-neighborhood of r* , hence d(r*) < 8. 

• If Nd{i"*) = 3 (implying f(r*) < 3), then, by Lemma [7] there are at 
most 16 petals in the 2-neighborhood of r*, hence d(r*) < 16. 

• If Nu(r*) = 2 (implying f(r*) < 2), then, by Lemma [BJ there are at 
most 20 petals in the 2-neighborhood of r*, hence d(r*) < 20. 

• If Nd{i"*) = 1 (implying f(r*) < 1), then, by Lemma [SJ there are at 
most 21 petals in the 2-neighborhood of r*, hence d(r*) < 21. 

For k = 1, ... ,4, let nk denote the number of relievers in G containing 
exactly k neighbors in D, so that \R'\ = Y^k=i n k- Let also C4 C C' be 
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the set of cores having at least one reliever with exactly 4 neighbors in D, 
and C 3 = \C \ C4I be the set of cores whose relievers have at most 3 vertices 
in D. Finally, we let £3 = |C 3 | and i 4 = |C^| , so that t = \C'\ = £3 + t±. 

Since G' is bipartite, the number of edges incident to C and the number 
of edges incident to R' are the same. Consequently, 

£ d(c*) = ^ d(r*), 

c*sC" r*eR' 

and 

2i 3 + 3i 4 <8n 4 + 16n 3 + 20n 2 + 21ni. (2) 

The same reasoning holds for the subgraph of G' induced by the cores 
in C 3 and their relievers in R', so 

2t 3 < 16n 3 + 20n 2 + 21n x . (3) 
By dividing both sides of ([3]) by 2 and adding it to ((21), we obtain 

63 

3(t 3 + *4) < 8n 4 + 24n 3 + 30n 2 + y n x . (4) 

As for the average p' of /(•) over the elements of G', we can write 

, 5(t 3 + i 4 ) + 4n 4 + 3n 3 + 2n 2 + n\ 
P = • 

h + *4 + "-4 + "-3 + "-2 + ".I 

Now, substituting (i 3 + 1 4 ) in the expression above by the upper bound 
obtained from (jK), we have 

. ¥n 4 + 43n 3 + 52n 2 + Wm 

P < • 

~~ ^n 4 + 9n 3 + lln 2 + 

Thus, p' is bounded by a multivariate rational function whose maximum 
is 43/9 = 4.777 . . ., achieved when n\ = n 2 = n 4 = 0. 

□ 

The following theorem summarizes the result from Section HI 

Theorem 13. Given the adjacency list representation of a unit disk graph 
with n vertices and m edges, we can find a 43 /9 '-approximation to the min- 
imum dominating set in 0(n 2 m) time. 
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Figure 6: Lower bound of 4.25 to the approximation factor of the modified 
graph-based algorithm. The coordinates of the centers are given in Table [5J 



(0,0), 

(±336577, 3647829), 

(±469471,882947), 

(±1043683,2940723), 

(±2657983,1789254), 

(±997564,4430244.0), 

(±5515150,2304292), 

(±2373785,4388387), 

(±515038,-857167), 



(0,45), 

(±3372414,3440722), 

(±2857376,5297889), 

(±2506389,2940723), 

(±1775036,1258725), 

(±4605648,790625), 

(±4605648,2787883), 

(±4657983, 1789254), 

(±999780,20942), 



(±3657983,1789254), 

(±3887452,5297889), 

(±892089,1789254), 

(±529919,5348048.0), 

(±5515150,1274216), 

(±1775036,2258725), 

(±3372414,4440722), 

(±4371043,4388387). 



Table 2: Coordinates of the centers of the disks in Figure El where all 
diameters are equal to 1000001. 



5 Conclusion and open problems 

We introduced novel efficient algorithms for approximating the minimum 
dominating set and minimum independent dominating set in unit disk graphs. 

On one hand, a linear-time algorithm was devised attaining a sub-5 
approximation factor, namely 44/9 < 4.889. Nevertheless, the best lower 
bound we know for the proposed algorithm is 4.8, which corresponds to the 
unit disk graph given in Figure Closing this gap would likely require 
the development of new tools to prove that certain graphs are not unit disk 
graphs, whereupon computer generated proofs may be useful. 

On the other hand, an enhanced approximation factor of 43/9 < 4.778 
was obtained by allowing for more local replacements, yet a lower bound of 
4.25, corresponding to the unit disk graph given in Figure [6] (with coordi- 
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nates in Table [2]), is the best we are aware of. Notwithstanding the 0(n 2 ) 
factor increase on its time complexity, such modified algorithm is still by 
far faster than, say, the state-of-the-art 4-approximation algorithm, which 
requires a geometric model as input and runs in 0(n 9 ) time. Moreover, since 
the number of (weak) reductions that are performed remains linear, it may 
be possible to conceive either a refined analysis or a smarter way to execute 
them. 
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